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^ ■ We analyze the Mass Varying Neutrino (MaVaN) scenario. We consider a minimal 

^ . model of massless Dirac fermions coupled to a scalar field, mainly in the framework of 

! finite temperature quantum field theory We demonstrate that the mass equation we 

\ find has non-trivial solutions only for special classes of potentials, and only within 

■ certain temperature intervals. We give most of our results for the Ratra-Peebles 

Dark Energy (DE) potential. The thermal (temporal) evolution of the model is 
PLh. analyzed. Following the time arrow, the stable, metastable and unstable phases 

Q ! are predicted. The model predicts that the present Universe is below its critical 

^ \ temperature and accelerates. At the critical point the Universe undergoes a first- 

^ I order phase transition from the (met a) stable oscillatory regime to the unstable rolling 

""^^ regime of the DE field. This conclusion agrees with the original idea of quintessence 

CT) . as a force making the Universe roll towards its true vacuum with zero A-term. The 

. present MaVaN scenario is free from the coincidence problem, since both the DE 

^ \ density and the neutrino mass are determined by the scale M of the potential. 

I Choosing M ^ 10~^ eV to match the present DE density, we can obtain the present 

neutrino mass in the range m ~ 10~^ — 1 eV and consistent estimates for other 
^ ■ parameters of the Universe. 

a\\ 
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I. INTRODUCTION 

Neutrino mass related questions are of great interest for particle physics as well as for 
cosmology (for reviews see Ref. [1] and references therein). Current upper limits on the 
sum of neutrino masses from cosmological observations are of the order of 1 eV [1-0] , while 
neutrino oscillations give a lower bound of roughly 0.01 eV [1,0], making neutrino mass an 
established element of particle physics. Furthermore, understanding the origin of neutrino 
mass opens a window into understanding physical processes beyond the standard model of 
particle physics 



It is now well established that about seventy four percent of the Universe is comprised 
of dark energy (DE) (for reviews see Ref. \)A\ and citation therein). The present stage of 
evolution of the Universe is governed by this dominant DE contribution, and the Universe 
experiences an accelerating expansion (l2L Il3l]. The nature of DE is still unknown, and it is 
one of the major questions of modern cosmology. There are, broadly speaking, three major 
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possibilities proposed to explain the DE fill]. Most straightforwardly, and in good agreement 
with the current observational data, it can be present just as the cosmological constant 
Secondly, the DE can be accommodated in some framework of the modified non-Einsteinian 
gravity theories (see, e.g., Refs. [3, [HI). And lastly, following the original proposals 17 



on the DE originating from a scalar field action similar to the infiaton field, there has been 
a lot of activity in constructing and analyzing various trial scalar field Lagrangians to model 
the DE [13]. Note, that it is even unclear what kind of scalar field potential governs the 
infiationary expansion of the Universe [18], and as the result, the eflFective quantum field 
that adequately describes infiation is still under debate jl9|. A similar observation can be 
drawn from analyzing many potentials proposed for the DE action [13[ |. 

On the other hand, several cosmological and astrophysical observations imply that about 
twenty two percent of the Universe consists of dark matter (DM) pjj], if we admit the general 
relativity theory of gravity. Most probably DM is formed through massive weakly interacting 
particles (WIMPs), and the nature of these particles is also still unknown. There are several 
recent observations performed by PAMELA and GLAST missions which indicate DM 
particle annihilations [i21i]. Recently it was proposed that both these observations could be 
used to test baryogenesis [22| which is one of the important problems of the standard particle 
physics model. 

Another puzzling question in modern cosmology is the coincidence problem - the density 
of DE is comparable to the present energy density of DM. In turn, the latter is comparable 
(within the order of magnitude), to the energy density of cosmological neutrinos [H, H])- Is 
there a mechanism explaining this coincidence? A very convincing answer to this question is 
given by the mechanism of DM mass generation via various types of DM-DE couplings, rang- 
ing from Yukawa to more exotic ones, fisl-lisjl The mass of the DM particle in this approach 
is naturally time-dependent, and they were coined Varying Mass Particles (VAMPs). Vari- 
ous DE-DM interaction models have been constrained by observations of Supernovae type 
la |29| . the age of the Universe js^-S], Cosmic Microwave Background (CMB) anisotropics 
|33|, l3J], and Large Scale Structure (LSS) formation |35l] . 



Pardon. Nelson and Weiner elaborated on the VAMP mechanism in the context of neu- 
trinos. [36]^. In their model the relic neutrinos, i.e., fermionic field(s), interact with a scalar 
field via the Yukawa coupling. If the decoupled neutrino field is initially massless, then the 
coupling generates a (varying) mass of neutrinos in this DE-neutrinos model. This mass 
varying neutrino (MaVaN) scenario is quite compelling, since it connects the origin of neu- 
trino mass to the DE, and solves the additional coincidence problem of why the neutrino 



mass and DE are of comparable scales [38|]. (For more on the coincidence, see, e.g. [39]). To 
consider neutrinos as particles which get their mass through the coupling is attractive for 
particle physics, as well as for its cosmological consequences. However there are significant 
issues that have to be resolved for the sake of viability of the MaVaN scenario. Most notably, 
it has been shown [4Q] that the model of Ref. [36] suffers from a strong instability due to 
the negative sound speed squared of the DE-neutrino fiuid (see also [41]). 

Any DM-DE coupling induces observable changes in large scale structure formation jiil ]. 
The main reason for this is due to the presence of additional DM contributions (perturba- 
tions) in the equation of motion which determines the dynamics of the scalar field. The 
changes in the dynamics are drastic when massive neutrinos are coupled to DE jiO]. In 
this case the squared sound speed of the DE-neutrino fiuid defined as = 5P/5p^ (where 



^ The DE-neutrino coupling and the baryogenesis constraints have been also studied also in Ref. [s^l 
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5 represents the variation, and P and p are pressure and energy density of the DE-neutrino 
fluid) is negative. The negative squared sound speed results in an exponential growth of 
scalar perturbations. [43>-46j 

After the critique in Ref. |40l] . t he issue of stability of the DE-neutrinos fluid has been 
addressed by many authors j4l|, l46l-l53| . Various physical assumptions were made in those 
references in order to avoid the exponential clustering of neutrinos. In particular, to achieve 
stability, proposals were put forward to make the DE-DM model more complicated, e.g., by 



extending it to a multi-component scalar fleld, or by promoting its supersymmetry. |49|, [51 
We however are not inclined to pursue this line of thought and will explore the simplest 
possible "minimal" model. As we will demonstrate, the occurrence of the instability in the 
coupled DE-neutrinos model is meaningful, and we will explore the physical implications of 
this phenomenon. Note that Wetterich and co-workers \0\ have already analyzed various 
implications of the instability in the MaVaN model on the dynamics of neutrino clustering. 

In this paper we re-address the analysis of the DE-neutrinos coupled model. What is 
really new in our results, to the best of our knowledge, apart from a consistent equation 
for the equilibrium condition, is the analysis of the thermal (i.e. temporal) evolution of the 
MaVaN model and prediction of its stable, metastable and unstable phases. The analysis 
of the dynamics in the unstable phase results in, for the flrst time in the framework of the 
MaVaN scenario, a picture of the present-time Universe totally consistent with observations. 
Our flndings are in line with the original proposal (l6,|l3] of the DE potential (quintessence) 
to model the Universe slowly rolling towards its true vacuum (A = 0). As it turns out, the 
present Universe, seen as a system of the coupled DE (quintessence) fleld and fermions 
(neutrinos) is below its critical temperature. It is similar to a supercooled liquid which has 
not crystallized yet: its high temperature (meta)stable phase became unstable, but the new 
low-temperature stable phase (A = 0) is still to be reached. The Afshordi-Zaldarriaga-Kohri 
instability corresponding to c| < is just telling us this. 

The rest of the paper is organized as follows: In Section [TTl we give the outlook of the 
model and formalism applied and derive the basic equations for the coupled model. In 
Section IIIII we present the qualitative analysis of the equation which yields the fermionic 
(neutrino) mass. Section HVl contains analysis of the coupled model with the Ratra- Peebles 
DE potential at equilibrium. The dynamics of the model applied to the whole Universe is 
studied in Section [Vl The results are summarized in the concluding Section IVll 



II. MODEL AND FORMALISM. BASIC EQUATIONS 

A. Outlook 

In this paper we focus on the case when the scalar fleld potential U{(f) does not have 
a non-trivial minimum, and the generation of the fermion mass is due to the breaking of 
chiral symmetry in the Dirac sector of the Lagrangian. A non-trivial solution of the fermionic 
mass equation is a result of the interplay between the scalar and fermionic contributions. 
We consider the most natural and intuitively plausible Yukawa coupling between the Dirac 
and the scalar flelds. 

The key assumption is that the fermionic mass generation can be obtained from mini- 
mization of the thermodynamic potential. That is, the coupled system of the scalar bosonic 
and fermionic flelds is at equilibrium, at least at some temperatures. This will be analyzed 
below more speciflcally. We assume the cosmological evolution, governed by the scale factor 
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a(t) to be slow enough that the coupled system is at equilibrium at a given temperature 



T(a). Then the methods of thermal quantum field theory [5J, l55|] can be applied. 

This problem is rather well studied with quantum field theory and statistical physics 
in different contexts [54-5^. The major conceptual difficulty in applying quantum field- 



theoretical methods for the dark-energy scalar field is the lack of "well-behaved" potentials 
interesting for cosmological applications. For instance, a class of the very popular inverse 
power law slow-rolling quintessence potentials fis'l are singular at the origin. Consequently, 
the field theory should be understood as a sort of effective theory, and we plan to address 
this issue more deeply in our future work. 

As far as the fermionic sector of the theory is concerned, one needs to distinguish two 
different cases pertinent for neutrino applications: 

(i) an equal number of fermions and antifermions, i.e., zero chemical potential /i = 0; 

(ii) a surplus of particles over antiparticles, and small non-zero chemicalpotential. 

For the bounds on the neutrino chemical potential, see Refs. l53|. If experiments 
confirm neutrinoless double beta decay, i.e., that neutrinos are Majorana fermions, then 
the lepton number is not conserved [8], and one cannot introduce a (non-zero) chemical 
potential. Then case (i) above is applicable, proviso that the Majorana fields are utilized 
instead of the Dirac ones. For the case (i) with Dirac fermions the ground state corresponds 
to a complete annihilation of fermion-antifermion pairs, i.e. the fermions completely vanish 
in the zero-temperature limit. 

Assumption of the fermion-antifermion asymmetry and (conserving) particle surplus, i.e., 
of a non-zero chemical potential, results in the fermionic contributions which survive the zero- 
temperature limit. However the smallness of the zero-temperature contribution renders this 
issue rather academic. Indeed, for the neutrinos we are interested in this study, by assuming 
the maximal particle surplus Uo ^ 115 cm~^, one g ets the Fe rmi momentum ^ 3-10~^ eV. 
For m ^ 10"^ eV, one obtains /i(r = 0) = = ^/l^^^^m^ = m + 0(10"^ eV). This results 
in a non-trivial vacuum with the particle surplus frozen within an extremely narrow Fermi 
shell m < e < Sf- Thus, trying to grasp the essential physics in this study from possibly 
the simplest "minimal model", we assume the fermions to be described by a Dirac spinor 
field with zero chemical potential. 

In this work we will use the standard methods of general relativity and finite-temperature 
quantum field theory extended for fields living in a spatially fiat Universe with the 
Friedmann-Lemaitre-Robertson-Walker (FLRW) metric where the line element is ds^ = 
dt^ — a^{t)d'x^ . Here t is the physical time and a{t) is the scale factor, which can be obtained 
from the Friedmann equations [1, l3 



H'{t) = - = ^Ptot , (1) 



a J 3 

Hit) + H\t) = - = -^(Aot + 3Ptot) . (2) 
a S 



Eqs. (in)-® also lead to the continuity equation 



Aot + — (Ptot + i^tot) = . (3) 



a 



Here the dot represents the physical time derivative and ptot and Ptot are the total energy 
density and pressure of the Universe. In accordance with the (standard) ACDM model, the 
Universe is assumed to consist of (1) DE, (2) cold DM (CDM) made of weakly interacting 



5 



massive particles, presumably Mdm > 1 ^ 10 GeV, (3) photons, and (4) baryons. The 
DM and baryon density parameters today are Qbm = PBuitnow) / Pa ^ 0.22 and = 
Pb{Uow)/pcr ^ 0.04. Here p^r = 3H^/{87rG) = 8.1h^ x 10"^^ GeV^ is the critical density 
today, tnow defines the current time, Hq = 2.1/ix 10""^^ GeV is the present Hubble parameter, 
G is the Newton constant, and h ^ 0.72 is the Hubble parameter in units of 100 km/sec/Mpc. 
The photon contribution to the energy density today can be neglected. The fiatness of the 
Universe leads to the relative energy density of the DE-neutrino coupled fiuid Q^jy ^ 0.74. 
To ensure the accelerated expansion of the Universe today, the r.h.s. of Eq. (Ej) must be 
positive at t = tnow 

In this paper we will not assume the existence of the cosmological constant A, as the 
ACDM model suggests. Instead we accept the hypothesis of the dynamical dark energy 
described by a scalar field. This is a bold assumption and a highly debatable issue. We 
vindicate our approach a posteriori by the consistent picture we arrive at the end. For a 
review and/or alternative approaches, see, e.g., Refs. fisl, l58|, l59| . The massless neutrinos 
are described by the conventional Dirac Lagrangian. The resulting model is given by the 
coupled Dirac and scalar fields. The grand thermodynamic potential of the coupled model 
can be derived from the euclidian functional integral representation of the grand partition 
function. The dynamics of the coupled model is governed by the Friedmann equations. 

Throughout the paper we use natural units where h = c = kB = I- 



B. Bosonic Scalar Field 



The bosonic scalar field Hamiltonian in the FLRW metric reads as jsl, [60| 
Hb = J a^d^x 



(4) 



where the comoving volume V = J d?x^ while the physical volume Vphys = o?{t)V . Since 
this field does not carry a conserved charge (number), the chemical potential /x = 0. The 
grand partition function in the functional integral representation: 



Vipe 



with the bosonic euclidian action 



(5) 



(6) 



where ip (/^(x, r). 

It is instructive to find the partition function of the free scalar field U{lp) = ^M^^p'^ 
following the methods explained by Kapusta and Gale [54^ for the case of the Minkowski 
metric. Rescaling of the field 

^ = a^/V (7) 

changes the partition function (IS]) by a thermodynamically irrelevant prefactor. The func- 
tional integration over (p of the Gaussian action gives 



log Zb = -V 



(2^ 



PJm^ + k^o? + log f 1 - e-P^/Ml+kVa^ 



(8) 
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Then the density (with respect to the physical volume) of the thermodynamic potential is 
given by 



= log Zb = -Pe 



where e = ^/M^~^^^ and Pb is the pressure due to the bosonic field. 



(9) 



C. Free Dirac Spinor Field 

The Dirac Hamiltonian in the FLRW metric is [60] 

Hd = J a^d^x ^( - ^7 • ^ + • (10) 
The grand partition function is given by the following Grassmann functional integral: 

Zn = Tre-^(^-^^) = J V^V^e'^^ (11) 

where the conserved charge (lepton number) operator Q = J a^cPx^p'^i/j and the euclidian 
action ^ 

SE = drj a(t)3d3x V;(x,r)(7°|--^7- V + m-^7°)^(x,r). (12) 

By rescaling the Grassmann fields ([71) and using the standard techniques [54], we get the 
thermodynamic potential density (pressure) as a function of the chemical potential and 
temperature: 

where 

£{k) = Vm2 + fc2 , (14) 

and £± = e{k) ± ji. The first term on the r.h.s. of Eq. ( IT3l) corresponds to the vacuum 
contribution to the thermodynamic potential (pressure): 

J^A^) (15) 
Introducing the notation for the Fermi distribution function 

Eq. ( fT3l) can be brought to the following form: 

-nn = Pn = Po + ^,l 7(^M^-) + M^+)] (17) 
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D. Coupled Model: Scalar Field and Dirac Massless Fermions 

Let us consider a scalar bosonic field interacting via a Yukawa coupling with massless 
Dirac fermions. The euclidian action of the model in the FLRW metric reads: 

S = SE + SE\^^, + gjyTj a^d^x (18) 

The path integral for the partition function of the coupled model is: 

Z = j V^Vi^V^e-^ (19) 

The Grassmann fields can be formally integrated out resulting in 

Z = J P(^e-^(^) = Jv^exp[-S§ + logDet^((^)] , (20) 

where the Dirac operator 

b{^) = 7 V - -7 ■ V + ^^(x, r) - (21) 
or a 

The thermodynamic potential Vt of the model f fTSl ) at tree level can be found by evaluating 
the path integral f l20l ) in the saddle-point approximation. Assuming the existence of a 
constant (x, r)-independent field (f)c which minimizes the action S{(^)^ the term logDetZ? 
can be evaluated exactly, and fermionic contribution to the thermodynamic potential is 
given by Eqs. (fT3l) or (fTTI) with the fermionic mass 

m = g(j),. (22) 

The bosonic contribution to the partition function in this approximation is simply Z oc 
exp[— /3a"^y{7((/)c)] ■ The thermodynamic potential density is given then by 

n{4),) = u{4>c) + ^D{4>c) . (23) 

Self-consistency of the employed saddle-point approximation naturally coincides with the 
condition of minimum of the thermodynamic potential at equilibrium (at fixed temperature 
and chemical potential): 

— ^ = , (24 

>P=4>c 

and 



> , (25) 



Note that a non-trivial solution (j)c of Eq. IHM (if it exists) is called the classical field: it 
is the average of the bosonic field, i.e., (j)c = {(f)- Eqs. f l22l[23l[24l ) can be brought to the 
equivalent form: 

U'{<P,)+gp, = , (26) 
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where the scalar fermionic density (a.k.a. the chiral density) Ps is given by the foUowing 
expression: 

Ps^^ = ^ = Po + -,l [nFie.) + n,ie^)] , (27) 

and N = J d^xipip. Here po stands for the vacuum contribution to the chiral condensate: 

p^^^ = _!iL (28) 

dm TT^ Jo ^(^) 

Note that even if the time, i.e., a(t), does not enter explicitly in the equations for the 
thermodynamic quantities of the coupled, fermionic or bosonic models f l9lll3ll23ll26l[27l) . and 
they look like their counterparts in a flat static Universe, such parameters as, e.g., the 
temperature and chemical potential in those equations are time-dependent, i.e., T = T{a) 
and ji = /x(a). The particular form of the dependencies T{a) and /x(a) must be determined 
from the Friedmann continuity equation ([3]) which relates the energy density p{T) and 
pressure P{T) to the evolution of a(t)[i, [13]. In addition, the fermionic mass m oc (pc'm the 
coupled model is also time varying, since the time enters into (pc (l26l) via T, /x, and all three 
functions m(a), T(a) and //(a) are governed by the Friedmann equations ([I]|2]|3]). 

The present theory works consistently for the physical quantities (bosonic or fermionic) 
measured with respect to their vacuum contributions. So, in the rest of the paper we 
will employ the thermodynamic quantities with subtracted vacuum contributions, keeping 
however, the same notations, e.g.: 

Qd^^d-^o ^Pd^ Pd- Po ^Ps^ Ps- Po ' (29) 

Then, according to Volovik [g]]], the pressure and energy of the pure and equilibrium vacuum 
is exactly zero. (The renormalization of the vacuum terms is, of course a very subtle issue. 
There are alternative approaches to this problem known from the literature. See, e.g.. 



MM-) 



III. ANALYSIS OF THE MASS (GAP) EQUATION: GENERAL PROPERTIES 

In cases interesting for cosmological applications, the scalar fleld potential U{(f) does not 
have a non-trivial minimum, and the generation of the fermion mass (i.e. a solution of IHM 
< (/)c < oc) is due to the interplay between the scalar and fermionic contributions to the 
total thermodynamic potential (l23l) . 

From now on we adapt our equations for the case of equal number of fermions and 
antifermions and /x = 0, as discussed in Sec. Ill A[ Keeping in mind the neutrinos, we assume 
an extra flavor index of fermions with the number of flavors 5. (For neutrinos 5 = 3.) We 
also assume the flavor degeneracy of the fermionic sector. 

Before proceeding further, we need to make some important observations regarding the 
behavior of the coupled model in two limiting cases. Assuming that a non-trivial solution 
of IHM with flnite m exists, the fermionic contribution to the thermodynamic potential 
(pressure) (fTTI) can be written as: 

2q 

-nr) = PD = , ^ = , (30) 
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where the integral defined as 



X,{k) = / fl'^' dz (31) 



can be evaluated analytically in two cases: 



W = (32) 
3K^K2{K) + 0{e-^^) , K>1 



where K^^^x) is the modified Bessel function of the second kind. 
In the (classical) low-temperature regime 



m 

/5m = - » 1 (33) 

the above equation results in 

-nn = Pr, = A^K^i^m) + 0(e-2/^'") . (34) 

To leading order 

-^d = Pd^ ^T{Tmfl\-^l^ . (35) 
The chiral condensate density ( 1271) 

can be also evaluated in the low-temperature limit as 

Ps = ^Kr{/3m) + 0{e-'^n , (37) 



which gives to leading order 



^ ^^Tmyi'e-^'' . (38) 

In this limit the fermions enter the regime of a classical ideal gas. Indeed, the fermionic 
particle (antiparticle) density 

s r - (/9m)2)l , 

in the low-temperature limit yields 

n± = ^i^2(/3m) + Oie-^^n , (40) 

and to leading order: 

^ -I—^Trnf^e-""/^ . (41) 



\/27r3/2 
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We see from Eqs. f l34ll40l) that up to terms 0{e ^Z^^), the fermions satisfy the ideal gas 
equation of state 

PD^{n^ + n_)T , (42) 
and the chiral density is equal to the total particle density n. 

pg ^ n = + . (43) 

In the (ultra-relativistic) high-temperature regime 

« 1 (44) 



one obtains 



T 



-nr, = P^^'I^T'-^{mTf . (45) 



To leading order the chiral condensate is 



while the particle density is 



Ps - ^mT^ , (46) 
6 



n± ^ . (47) 



Now we can make some general observations of the fermionic mass generation in the 
coupled model: 

(i) It is obvious from the sign of Ps (cf. [2711361) that non-trivial solutions of l[2E\) are 
impossible for a monotonically increasing potential U{(f). That rules out some popular 
potentials, e.g., U (x log(l + (f/M) [H, 3y| for this Yukawa- coupling driven scenario of the 



mass generation. 

fii) The monotonously decreasing slow-rolling DE potentials ([l6|, [l3| and for reviews, 
see j3, EBtj), e.g., U oc Lp~^ oi U oc exp[— A^^^], do have a window of parameters wherein 
non-trivial solutions of ( l26l) exist. As we can see from ( l38l) . for those decreasing potentials 
the mass equation (l26l) always has a trivial solution m = g(j)c = oc for the minimum of 
the thermodynamic potential ( l23l) . ^ This solution corresponds to a "doomsday" vacuum 
state [61], when the Universe reached its true ground state with zero dark energy density 
and completely frozen out fermions. A non-trivial solution of ( l26l) . corresponding to another 
minimum of the potential ( l23l) . is totally due to the fermionic contribution. Since the latter 
freezes out in the limit T ^ 0, it is clear qualitatively that such a solution < m < cx) can 
exist only above a certain temperature. For a more quantitative account of these phenomena 
we need to assume some specific form of the DE potential. This will be done in the following 
section. 

(Hi) To explain the differences between the present study and earlier related work on 



mass varying fermions (see [23,, l2J, |36|, l40|| and more references there) , some clarifications 
are warranted. It is usually assumed in the literature that the low-temperature regime 
formulas are applicable, and according to ( ll3l) Ps = n. The approximation for ( l26l) then can 
be written as dU /dm + n = 0. The latter is interpreted as a result of minimization of some 



^ Recall that the grand thermodynamical potential is equal to the free energy for the case ja = 0. 
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effective potential [/gfr = U + nm with fixed n, which always has a non-trivial minimum 
< m < oc for the class of decreasing potentials [/, see, e.g., |23|, l24l- It turns out that such 
an approximation changes the picture qualitatively. 

In what follows, we explore in detail the predictions of the consistent mass equation ( l26l) 
on the mass varying scenario for the coupled model with a specific DE potential ansatz. 



IV. COUPLED MODEL WITH THE RATRA-PEEBLES QUINTESSENCE 

POTENTIAL 

A. Mass Equation and Critical Temperature 

Now we analyze in detail our coupled model for a particular choice of U{lp)^ the so-called 
Ratra-Peebles quintessence potential |1( 



= , (48) 

where a > 0. It is convenient to introduce the dimensionless parameters 

M qm ^ VL 

^ = ^^T^^^^W ^^^^ 

Then the mass equation f l26l ) can be written as: 

g"A''+^=XM, (50) 



25 

where we introduced 



= K'^^' I ^^^^dz . (51) 



According to the relation Eq. ( I22l) between the fermionic mass m and the classical field, we 
get m = Tkc^ where Kc is the solution of Eq. ( I50l) corresponding to the minimum of the 
thermodynamic potential which reads now as (cf. Eq. (IHTI)): 

"- = «°(7)°-j5ilW. (52) 

The dimensionless Yukawa coupling constant g ^ 1. To reduce the number of model pa- 
rameters we can set 5^ = 1. This is equivalent to the simultaneous rescaling gip (f and 
Mg^-^ M. ^ For simplicity, we also restrict the number of fiavors 5=1. 
We define the mass of the scalar field as: 



(53) 



^ One can check this scaling also holds for the dynamics of the model, considered in Section|Vl In particular, 
the neutrino masses do not depend on the value of ^. To avoid cluttering of notations we will drop tildes 
in the rescaled parameters. 
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In terms of the dimensionless parameters it reads 

i=^^(f)^ ,54) 

It is important to realize that the integral X(y,{n) on the r.h.s. of the mass equation is 
bounded. The quantitative parameters of the function depend on a, but its shape 

is always similar to the curve shown in Fig. [T]for a = 1. So, there exists a maximal Acrit 
(critical temperature Tcrit) such that for A > Acrit {T < Tcrit) only a trivial solution m = oo 
exists, and the stable vacuum has zero energy and pressure. 




FIG. 1: (Color online) Left: Graphical solutions of the mass equation (l50|l for different values 
of A = M/T (a = 1). Right: dimensionless density of the thermodynamic potential (l52]) . The 
thermodynamically stable solutions of Eq. (|50|) indicated by the large dots correspond to the minima 
of the potential. The arrows indicate the unstable solutions of the mass equation, corresponding 
to the maxima of the potential. 

The mass equation Eq. f l50l) is solved numerically for various values of its parameters, and 
the characteristic results are shown in Fig. [H The numerical results can be complemented 
by an approximate analytical treatment of the problem. The latter turns out to be quite 
accurate and greatly helps in gaining intuitive understanding of the results. 

It is easy to evaluate Ia{i^) to leading order: 

(55) 
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For the critical point where X^^^crit) = 0, we obtain: 

5 , , 

/^crit ^ ^ , Ty = a + -] (56) 

lai^crit) ^ yi^'e"'- (57) 

The most important conclusion we draw from Fig. [T]is that there are three phases in the 
modeVs phase diagram. We analyze each of them in the following subsections. 



1. Stable (massive) phase: A < Ao (To < T < oc) 

In this range of parameters the equation ( l50l ) has two nontrivial solutions. The root 
hie < indicated with a large dot in Fig. [1] (case a) gives the fermionic mass and corresponds 
to a global minimum of the potential. So it is a thermodynamically stable state. In this 
phase ^{tvc) < 0, so the pressure is positive P > 0. Another non-trivial root of ( l50l) 
corresponds to a thermodynamically unstable state (maximum of Q indicated with an arrow 
in Fig. d]). There is a trivial third root of the mass equation = oc. At these temperatures 
it corresponds to the metastable vacuum state Q = 0. 

In the high-temperature region of this phase where A <C 1 the fermionic mass is small 
(see Fig.Ej): 

^«(yfef)^ocr-A (58) 

The fermionic contribution to the thermodynamic potential is dominant, and it behaves to 
leading order as the potential of the ultra-relativistic fermion gas (cf. Eq. (Il5l)): 

n = -P = -^T^ + O(T^) . (59) 
180 

One can check that the subleading term in the above expression combines the DE potential 
contribution and the first fermionic mass correction, which are both of the same order. 

It is important to stress that in this coupled model with the slow-rolling potential Eq. ( Il8l) . 
the mass generation does not follow a conventional Landau thermal phase transition scenario. 
There is no critical temperature below which the chiral symmetry is spontaneously broken 
and the mass is generated. Instead the mass grows smoothly as /^c oc A«+2, albeit starting 
from the "point" T = oc. From physical grounds we expect the applicability of the model 
to have the upper temperature bound: 

T < Trd , (60) 

where Trd is roughly the temperature of the boundary between infiation and the radiation- 
dominated era. The high-temperature result ( l59l ) shows that the stable massive phase of 
the present model can indeed be extended up to those temperatures. 

The scalar field and fermionic masses demonstrate opposite temperature dependencies. 
The scalar field is "heavy" at high temperatures: 

~ V^y^^ , A < 1 , (61) 

however its mass decreases together with the temperature. In contrary, the fermionic mass 
m monotonously increases with decreasing temperature. The exact numerical results for the 
two masses are shown in Fig. [2] 
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FIG. 2: (Color online) Masses of the fermionic and scalar fields (m and m^^ resp.) as functions of 
A = M/T, a = 1. At A > Acrit {T < Tent) the stable phase corresponds to m = oc and m^^ = 

2. Metastable (massive) phase: Ao < A < Acrit (Tcrit < T < To) 



Upon increasing A we reach a certain value Ao corresponding to a critical temperature To 
when the thermodynamic potential has two degenerate minima Q{hio) = P{i^o) = f^(oc) = 0. 
This is shown in Fig. [1] (case b). After this point, when the temperature decreases further in 
the range Ao < A < Acrit (here Acrit stands for the maximal value of A when a non-trivial 
solution of the gap equation ( l50l) exists, see Fig. [1]), the two minima of the thermodynamic 
potential exchange their roles. The root k^c now becomes a metastable state with f^(/^c) > 0, 
i.e., with the negative pressure T(/^c) < 0, while the stable state of the system corresponds 
to the true stable vacuum of the Universe (giI] Q{oo) = P(oo) = 0. See Fig. [1] (case c). 
The system's state in the local minimum Q{tvc) is analogous to a metastable supercooled 
liquid. We disregard the exponentially small probability of tunneling of the fermions from 
the metastable state Q{tvc) into the vacuum state Q{oo) = o[l8|. Accordingly, the fermionic 
mass in this phase is determined by the root of ( l50l) . 

In the metastable phase /^c ^ 1, so by using Eqs. ( I52ll32ll50l) we obtain the potential: 



n 



R 



a 

Kip 



3q; ^ 
2^/ ■ 



(62) 



From the above result we can find the met ast ability point f](/^o) = as 



a f I 6 

^o^-(i + a/i + - 

a 



Expanding Xc,{n) near its maximum and using Eqs. 
Eq. f ISOl) . we obtain the following equation: 

(^c ~ ^crit)^ 



(63) 

along with the gap equation 



A \ 



2iy V Acrit > 

On finds from the above equation, e.g., how the mass approaches its critical value: 



merit - m oc 



^ crit 



1 



1/2 



(64) 



(65) 
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or the ratios of temperatures and masses at the metastable and critical points. These latter 
parameters are given in Table [B 

TABLE I: Masses, critical temperatures and potentials for various values of a. All the parameters 
used in this table are defined in the text. 



a 


^crit 

To 


^crit 


mo 


m-crit 

M 


<p 

M 


M 


^crit 


Pcrit 




1 


0.90 


0.91 


0.558 


3.86 


0.187 


1.10 


0.15 


0.84 


-0.18 


2 


0.95 


1.04 


0.70 


4.35 


0.130 


0.97 


0.02 


0.25 


-0.09 


4 


0.98 


1.44 


0.81 


4.52 


0.048 


0.70 


6 • 10-4 


0.02 


-0.03 


10 


0.99 


3.00 


0.91 


4.16 


2 • 10-3 


0.33 


7 • 10-s 


9 • 10-*5 


-0.008 



3. Critical point: A = Acrit {T — Tcrit ) and phase transition 

The critical point of the model corresponds to the case when the two roots of the mass 
equation Eq. (I50l) merge, and the minimum of the potential disappears. One can check that 
instead of the minimum this is an inflection point of the the potential, i.e., f^'^(/^crit) = 0. 
This situation is shown in Fig. [1] (case d). At this point the system is in the unstable state 
with the fermionic mass 

^crit 

TP, — = i^cnt ~ ■ (66) 

^ crit 

In particular, this implies that the fermions are non-relativistic at the critical temperature. 
From Eqs. (I57II50I) we flnd the critical parameter (see Table [I] for its numerical values) 



which allows us to evaluate the critical temperature 



M 

Tcrit = -r- . (68) 

^crit 



5 /Acrit \-^, .4 



We can also find the potential at Tcrit: 

^^cnt^^(^) (69) 

Thus, from the viewpoint of equilibrium thermodynamics at T = Tcrit the model must 
undergo Si first-order (discontinuous) phase transition and reach its third thermodynamically 
stable (atT < T^rit) phase corresponding to the vacuum Vt{n = oo) = P{k = oo) = 0. During 
this transition the fermionic mass given at the critical point by Eq. ( l66l) and the scalar field 
mass 



m 



~ V«(a + 1)( ^) ""'m (70) 



both jump to their values in the vacuum state m = oo and m^^ = 0. See Fig. [2l 
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However, the above arguments are based on the minimization of the thermodynamic 
potential (i.e. maximization of entropy) at equihbrium. To address the question of how 
such a system behaves as the Universe evolves towards the new equilibrium vacuum state, 
we need to analyze the dynamics of this phase transition. More qualitatively, we need to 
study how the particle at the point /^crit at the critical temperature (see Fig. [1]) rolls down 
towards its equilibrium at infinity. This issue will be addressed in Section [Vl 



B. Equation of State 

We define the equation of state in the standard form: 

P = wp, (71) 

where the total pressure in this model is obtained from Eq. ( I52l) . while the total energy 
density (p) and its dimensionless counterpart {pr) are determined by the following equation: 

Here we define the integral 



XM) ^ / / , (73) 



+ 1 

which can be evaluated in two limits of our interest: 

Uk) = (74) 
[ 3n^K2{n) + n^Ki{n) + C(e-2-) , n>l 

In the high-temperature region of the stable massive phase where A <C 1, the fermionic 
contribution is dominant, and the energy density to leading order is that of the ultra- 
relativistic fermion gas (cf. Eq. fjS^ ) 

P=^-^T' + 0{T^^) . (75) 

Thus, in this regime the model follows approximately the equation of state of a relativistic 
gas with w ^ ^. 

In the region /^c ^ 1 which includes the metastable phase and the critical point, we obtain 

'"'(s:) ^+°+^+2;3}. (76) 

and 

2 a_ 3q; 

W ~ '"'q^ (77) 

The last equation follows very closely the results of the exact numerical calculations shown 
in Fig. [3l At the critical point we evaluate 



17 




FIG. 3: (Color online) w = P/p for several values of a. At A > Acrit(c^), i-e., T < rcrit(<^) the 
equilibrium value w = —1 exactly. 



p„.^{^y{l + c.+ '-^}M', (78) 



and making a rough estimate, we get a lower bound: 

.5 I >_1 

2p{l + a + ^) - 4 



w ^ -- , , > --, Va>l. (79) 



Thus for any power law a > 1, the parameter w of this model at equilibrium cannot cross 
the bound w < — |, necessary for accelerating expansion of the Universe a > 0. 

At r < Tcrit we obtain the equilibrium value of w in the stable vacuum state from 
Eqs. 

w = lim = -1 . (80) 

So the true vacuum in this model corresponds to the Universe with a cosmological constant 
in the limit A ^ 0. 



C. Speed of Sound 

We define the sound velocity as 



dP dP^ 

-l = t = t^ (81) 

dt dA 



The relation (|79]) w{Tcrit) ^ — ^ holds for the model which contains only the DE-neutrino coupled fluid. 
In a more realistic model for the Universe, baryons and DM also contribute to the total energy density, 
and as a consequence (Tcrit) increases, see Sec. |Vl 
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where to obtain the second expression we used the fact that the time enters our formulas 
only through the temperature T(a(t)), so 



d dA d 
di ~ ~dtdA ' 

Let us first consider the temperatures T > Tc^t, i.e., A < Acnt- Then 

dp dp dp dn 
dK^dA^di^'dA 

where hi is related to A through the gap equation (l50l) : 

a + 4:Ia{i^c) 0^ + 4: fdlogIa{i^c^ ^ -1 



di^ 
dA 



— I^r — 



Note that for the pressure the following relation 

dP dP 



dtv 



holds, since 



dA OA 



dP 

Ok 



is just another form of the gap equation f l24l ). Thus 



dP 

dA 



dp_ I 



(82) 



(83) 



(84) 



(85) 



(86) 



(87) 



In the high-temperature regime A <C 1 (/^c <^ 1), it is even easier to use the explicit 
asymptotic expansions for P{A) and p(A) in the definition ( ISTI) instead of the above formula 
(1871) . A straightforward calculation gives the result 



1 , ^ 2(a+4) 

d ^ 6A -+2 , b> 

s 2 



(88) 



consistent with the earlier observation that for A <C 1 the model behaves as an ultra- 
relativist ic Fermi gas. 



In the case /^c ^ 1 we find 



Of + 4 /^c 
A u — Kr ^ 



and 



(89) 



(90) 



^ a(a + 4)(l + ^)- 

Everywhere at T > Tcrit, including the stable and metastable massive phases > 0, so the 
model is stable with respect to the density fiuctuations. The sound velocity vanishes in the 
limit T T^-^^ as 



Cs (X \/v — Kc . 

Qualitatively, the vanishing speed of sound is due to divergent 
point. 



(91) 

at the critical 
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FIG. 4: (Color online) The square of the sound velocity for several values of a. At A > Acrit(<^), 
i.e., T < Tcrit(<^) the equilibrium value = — 1 exactly. 

The above analytical results are in excellent agreement with the numerical calculations 
of cl from the formula (l87l) shown in Fig. [H At the temperatures T < Tcrit there is no gap 
equation relating tv and A, so the sound velocity is easily calculated to yield the value in 
the equilibrium vacuum state: 

dP 

cl = lim = -1 . (92) 

OA 

That what is expected for a barotropic perfect liquid with a constant where cl = w. 



V. DYNAMICS OF THE COUPLED MODEL AND OBSERVABLE UNIVERSE 

A. Scales and Observable Universe 

In order to make a connection between the above model results and the observable Uni- 
verse, we need to first conclude where we are now with respect to the critical temperatures 
To and Tent- As one can see from Table [J for a ^ 1, the model has To ^ TcHt ^ M. We 
identify the current equilibrium temperature of the Universe with the cosmic background 
radiation temperature T = 2.275 K = 2.4 • 10~^ eV. Then we see right away that we cannot 
be above the critical temperature of the coupled model, since: 

(i) assumption T > Tcrit leads to M < 10~^ eV, which in turn implies too small densities 
p ^ ^ 10~^^ eV^, i.e, four orders of magnitude less than the observable density; 

(ii) At T > Tcrit the equation of state has w > —j (see Fig. [3j), which is not even enough to 
get a positive acceleration a > 0, while the observable value w ^ —1. [HI 

So, the first qualitative conclusion is that we are currently below the critical temperature. 
The Universe has already passed the stable and metastable phases and is now unstable, i.e. 
it is in the transition toward the stable "doomsday" vacuum m = oo and = 0. 

Since at the temperature of metastability To = 0, the transition occurs somewhere be- 
tween the beginning of the matter-dominated era {Tmd ^ 16500 K ^ 1A2 eV) and now, 
i.e., 1.4 eV > Tcrit > T^ow ^ 2.4 • 10"^ eV. Because of Eq.(|68j) this inequality sives us the 
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possible range of the model's single parameter M: 

2.4.10"^eV <M< 1.4 eV. (93) 

As we will show in the following, other consistency checks of the model bring the upper 
bound of M much lower. 

B. Universe Before the Phase Transition 

In order to apply the results of the coupled model for the calculation of the parameters 
of the observable Universe, we need to incorporate the matter (we will just add up the dark 
and conventional baryonic matter together) and the radiation. Assuming a spatially flat 
Universe, the total energy density is critical, so 

Ptot P7,now/tt^ + pM,now/ci^ + P^u{^) Per , (94) 

where from now on we denote p^j^ the energy density of the coupled model given by Eq. ( f72l) . 
To relate our model's parameters to the standard cosmological notations, we assume that 
the temperature is evolving as that of the blackbody radiation, i.e., T = T^o^/a. Then 

M Ma M 
A = — = = . (95) 

-'now -'nowv-'- ' ^) 



We know that 



P7 = J^T' , (96) 

and we set the current density of the coupled scalar fleld to the observable value of the dark 
energy, i.e., 3/4 of the critical density: 

P^.,now = i ■ ^ ^ 31 ■ (10-3 eV)^ , (97) 

and 

PM,„ow ^ 4 ■ ^ ■ (98) 
The equations above allow us to plot the relative energy densities 

0# = p#/Aot (99) 

as functions of redshift (or temperature) up to the critical point, see Fig. [5l ^ In the 
high-temperature limit, the matter term is sub- leading and 

7^2 7 

Act ~ + P^. ~ - (1 + ^)T^ . (100) 



^ We apologize for some abuse of notations, but using the same Greek letter for the grand thermodynamic 
potential and relative densities seems to be standard now. Since these quantities are mainly discussed in 
different sections of the paper, we hope the reader will not be confused. 
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FIG. 5: (Color online) Relative energy densities plotted up to the current redshift (temperature, 
upper axis): Q^i^ - coupled DE and neutrino contribution; fi^ - radiation; ~ combined baryonic 
and dark matters. Parameter M = 2.39 • 10~^ eV {a = 0.01), chosen to fit the current densities, 
determines the critical point of the phase transition Zcr ^ 3.67. The crossover redshift z* ^ 0.83 
corresponds to the point where the Universe starts its accelerating expansion. 



In this limit, then 

n^r. = ^^ 0.636 , n^ = ^^ 0.363 , (101) 

which agrees well with the numerical results displayed in Fig. [5l At the critical point the 
matter strongly dominates and pu/ Pj,(^iy ^ 10^. 

The equation of state parameter of the entire Universe, Wtot^ is given by Ptot = '^totPtot- 
Since the matter contribution Pm = 0, then Ptot = + Pipu-> where P^ = ^p^ and the 
pressure of the coupled model P^jy is obtained from Eq. f l52l) . The numerical results of Wtot 
are given in Fig. [6l 



T (eV) 




-0.8 -\ \ \ \ \ \ \ \ \ 1 ' ' I -0.8 
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FIG. 6: (Color online) Equation of state parameter wtot = Ptot/ptot ioi M = 2.39 • 10 ^ eV 
{a = 0.01) plotted up to the current redshift (temperature, upper axis). 
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To analyze the dynamics of the coupled model we need, in principle, to go beyond the 
saddle-point approximation applied in the previous sections and solve the equation of motion: 

OQ 

^ + 3H^ + — = . (102) 

dip 

Above the transition point (T > Tcrit) the dynamics is quite simple. Let us analyze pertur- 
bations to the saddle-point solution of (Ei 



^(t) = (l>, + tp{t) . (103) 
Taylor-expanding the thermodynamic potential of the coupled model 

^=a;V + V(</>,)^2^ ... (104) 
dip 2 

with cj^ = f]'''((/)c), we obtain from ( 11021) the equation of a damped harmonic oscillator to 
the leading order: 

^ + 3H^ + uj^^ = . (105) 

So, the quintessence field Lp{t) oscillates around its saddle-point value (/)c with i/j{t) oc 
^lut-^m ^ The damping is very small, since as one can check 

> . (106) 

The violation of the above condition and breaking down of the oscillating regime occurs in 
the vicinity of the critical point, which is the infiection point of the potential (cj = 0). This 
is the well-known phenomenon of the critical slowing down near phase transition. Retaining 
the first non- vanishing term in ( 11041) . the equation of motion in the vicinity of the critical 
point reads: 

^ + + ^-^'\(\)^^^ = . (107) 

Neglecting the small damping term in this equation, its solution can be found analytically 
via a hypergeometric function. Since the explicit form of this solution is not very interesting 
at this point, we just emphasize the qualitative conclusion of the analysis: the fiuctuation 
'0(t) oscillates near the classical field in the stable (metastable) phase at T > T^it, and 
it enters the run- away (power-law) regime when T T^^'^^^\ 



C. Late-Time Acceleration of the Universe. Towards the End of Times 

The equilibrium methods are not applicable below the phase transition, and we study 
the dynamics of the model from the equation of motion (11021) together with the Friedmann 
equations ([I]l2|l3j) . Solution of the Dirac equations yields oc for the chiral density 
so the equation of motion (11021) at a < Ocrit reads: 

¥3 + 3//0 = -|^-p.,cHt(^)'. (108) 
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From the results of the previous section we evaluate the chiral density at the critical point: 

Ps,crit^a(^y^'M' . (109) 



The system of the integro-differential equations ( ll08l[Tl[2ll3l) was solved numericaUy. AU the 
quantities entering those equations are defined in the previous subsection, except that one 
needs to include the extra term ^Lp'^ in the computation of both ptot and Ptot- However 
the numerical results show that in the regimes of the parameters we are interested, the 
kinetic term can be safely neglected. Since the critical point of the model lies in the matter- 
dominated regime (cf. Fig. Ej), we start with the Hubble parameter H = 2/3t {a (x t^/^). 
At the latest times {z < 1) the Hubble parameter was determined self-consistently from the 
numerical solution of the Friedmann equations. 

We find numerically that the quintessence field (f{t) from the critical point to the present 
time oscillates quickly (with the period r ^ 10~^^ Gyr) around the smooth ("mean value"' 
solution (^(t), where the "mean" (f nullifies the r.h.s. of the equation of motion 



Relating the mean values with the physically relevant observable quantities, we can easily 
obtain the key results analytically. (They are checked against direct numerical calculations 
and found to be accurate within 5 % at most). Thus we get 

^ = ¥^crit • -T"^ , (110) 



l + Z 



3a 
a + 1 



= Pv'.crit • (t-T—- ) ' (111) 

where (pait ~ -^^M and p^,crit ~ (^^)"M^. Having a free model parameter M, we'll set 
it by matching the current density of the scalar field p^^now to the observable value of the 
DE density (l97jl. so 

M = (i/V^,now) ^ A-«T„"o?^ . (112) 

The exponent of the quintessence potential a is now the only parameter which can be varied. 
We define the time-dependent mass via the solution of the motion equation as m{t) = (^(t), 
thus obtaining an estimate for the present-time neutrino mass. Results for various a are 
given in Table [III There we also calculate the critical points parameterized by the redshifts 
Zcrit and the crossover points z*. The latter is defined as the redshift at which the Universe 
starts its late-time acceleration, i.e., where Wtot = the present time we find 

<o7 ^ -I ■ (113) 

As we infer from the data of Table [TTl the range of exponents a <C 1 corresponds to more 
realistic predictions for the neutrino mass [H, 0, [sl and for the crossover redshift z* |65| . 
For a = 0.01 we plot the evolution of the relative energy densities, the equation of state 
parameter, and the neutrino mass in Figs. [5|l6|I71 

We consider the quite artificial case of small quintessence exponent a as an ansatz crossing 
over smoothly from physically plausible potentials with, say, a = 1 or 2 to the logarithmic 
potential 

[/((^) = M^(l + alog^) . (114) 
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TABLE II: Model's parameters and observables for various a. All the entries in this table are 
defined in the text. 



a 


M (eV) 


rwnow (eV) 


^crit 




2 


9.75 • 10-2 


167 


392 


4.9 


1 


1.69 • 10-2 


44.6 


76.6 


2.3 


1/2 


6.33 • 10-3 


17.0 


27.7 


1.5 


10-1 


2.81 • 10-3 


2.82 


8.73 


0.93 


10-2 


2.39 • 10-3 


0.27 


3.67 


0.83 


10-3 


2.36 • 10-3 


0.027 


1.60 


0.82 




z+1 



FIG. 7: Neutrino mass m for M = 2.39 • 10"^ eV [a = 0.01) plotted up to the current redshift 
(temperature, upper axis). 



The latter often appears in various contexts 13|, 



VI. CONCLUSIONS 

In this paper we analyzed the MaVaN scenario in a framework of a simple "minimal" 
model with only one species of the (initially) massless Dirac fermions coupled to the scalar 
quintessence field. By using the methods of thermal quantum field theory we derived for 
the first time (in the context of the MaVaN or, even more broadly, the VAMP models) a 



^ The numerical results for small parameter a, as e.g. a = 0.01 taken for the plots, are virtually in- 
distinguishable for the cases of the Ratra-Peebles (|48|) or logarithmic (jll4p potentials. However the 
Ratra-Peebles potential at more "natural" a = 1,2 allows to probe the coupled fermionic-quintessence 
models in the search of heavy DM particle candidates. 
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consistent equation for fermionic mass generation in the coupled model. 

We demonstrated that the mass equation has non-trivial solutions only for special classes 
of potentials and only within certain temperature intervals. It appears that these results 
have not been reported in the literature on VAMPs before now. 

We gave most of the results for the particular choice of a trial DE potential - the Ratra- 
Peebles quintessence potential. This potential has all the necessary properties we needed for 
our task: it is simple, it satisfies the criteria we found for non-trivial solutions of the mass 
equation to exist, and it has only one dimensionfuU parameter- the energy scale M to tune. 
Also, at small values of the exponent a it effectively crosses over to the case of a logarithmic 
potential. We have checked that other potentials, e.g., exponential, lead to a qualitatively 
similar picture, but they have at least one more energy scale to handle, which we consider 
as an unnecessary complication at this point. 

We analyzed the thermal (i.e. temporal) evolution of the model, following the time 
arrow. Contrary to what one might expect from analogies with other contexts, like, e.g., 
condensed matter, the model does not generate the mass via a conventional spontaneous 
symmetry breaking below a certain temperature. Instead it has a non-trivial solution for 
the fermionic mass evolving "smoothly" from zero at the "point" T — oo. The scalar field 
is infinitely heavy at the same point. More realistically, we assumed the model is applicable 
starting at the temperatures somewhere in the beginning of the radiation-dominated era. 
We found that the DE contribution in this regime is subleading, and the model behaves as 
an ultra-relativistic Fermi gas at those temperatures. 

This regime corresponds to a stable phase of the model given by a global minimum of the 
thermodynamic potential Vt[ip). The temperature/time dependent minimum {(p) generates 
the varying fermionic mass m oc {^p). 

With increase in time, as the temperature decreases, the model reaches the point of 
metastability where its pressure (P) vanishes. From our estimates of the model's scales, 
we showed that this happens during the matter-dominated era of the Universe. At this 
point the system's ground state becomes doubly degenerate, and the potential 11 = at the 
non-trivial (finite) minimum {^p) as well as at the trivial vacuum (/^ = oc. 

Further on, at lower temperatures the system stays in the metastable (supercooled) state 
until it reaches the critical point where the local minimum of the thermodynamic potential 
disappears and it becomes an inflexion point. At this critical temperature the model un- 
dergoes a first-order (discontinuous) phase transition. At the critical point the equilibrium 
values of the fermionic and the scalar field masses discontinuously jump to the 'doomsday" 
vacuum state values m — oo and = 0, respectively. The square of the sound velocity and 
equation of state parameter w have the equilibrium values corresponding to the de Sitter 
Universe with a cosmological constant, i.e. = w = —1. It is worth pointing out that 

> in both the stable and metastable phases, and the sound velocity vanishes reaching 
the critical temperature from above. 

Since the equilibrium approach is not applicable below the critical temperature, we find 
parameters of the model from direct numerical solution of the equation of motion and the 
Friedmann equations. The single scale M of the quintessence potential is chosen to match 
the present DE density, then other parameters of the Universe are determined. We obtain a 
consistent picture: the phase transition has occurred rather recently at Zcru ^ 5 during the 
matter-dominated era, and the Universe is now being driven towards the stable vacuum with 
zero A-term. The expansion of the Universe accelerates starting from z* ^ 0.83. Setting 
a = 0.01 for M ^ 2.4 • 10~^ eV, we end up with the neutrino mass m ^ 0.27 eV. 
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The present results allow us to propose a completely new viewpoint not only on the 
MaVaN, but on the quintessence scenario for the Universe as well. The common concerns 
about the slow-rolling mechanism for the DE relaxation toward the A = vacuum are related 
to the question of what is the mechanism to set the initial value of the scalar field Lp where 
it evolves (rolls down) from. Our results demonstrate that up to recent times (i.e. above the 
critical temperature) the quintessence field was locked around its average (classical) value 
(99). Its value is determined by the scale M and the temperature. The average {^p) gives 
the fermionic mass at the same time. The scalar field is rigid (i.e. massive), although it 
softens (i.e., its mass decreases) as the system approaches the critical temperature. Above 
the critical temperature the scalar field can only oscillate around its equilibrium value (99). 
At the critical point the minimum of the thermodynamic potential becomes the infiexion 
point, the scalar field looses its rigidity (mass). Then the field can only roll down towards 
the new stable ground state = at (/p = oc. So physically, the critical point corresponds 
to the transition of the Universe from the stable oscillatory to the unstable rolling regime. 

A more sophisticated numerical study of the kinetics after the critical point is warranted 
in order to address such issues as the detailed description of the crossover between diflFerent 
regimes, and the clustering of neutrinos. These and some other questions are relegated to 
our future work. 

Acknowledgments 

We highly appreciate useful comments and discussions with D. Marfatia, B. Ratra, and N. 
Weiner. We are grateful to N. Arhipova, D. Boyanovsky, R. Brandenberger, O. Chkvorets, 
H. Feldman, A. Gruzinov, L. Kisslinger, the late L. Kofman, S. Lukyanov, and U. Wi- 
choski for helpful discussions and communications. We thank the anonymous referee for 
constructive criticism and comments which stimulated us to undertake deeper analyses of 
the model, and especially of its dynamics. G.Y.C. thanks the Center for Cosmology and Par- 
ticle Physics at New York University for hospitality. We acknowledge financial support from 
the Natural Science and Engineering Research Council of Canada (NSERC), the Laurentian 
University Research Fund (LURF), Scientific Co-operation Programme between Eastern 
Europe and Switzerland (SCOPES), the Georgian National Science Foundation grants # 
ST08/4-422. T.K. acknowledges the support from NASA Astrophysics Theory Program 
grant NNX10AC85G and the ICTP associate membership program. A.N. thanks the Bruce 
and Astrid Mc Williams Center for Cosmology for financial support. 



[1] A. D. Dolgov, Phys. Atom. Nucl. 71, 2152 (2008) |arXiv:0803.388 7 [hep-ph]]; Physics Reports 
370, 333 (2002). 

[2] J. Lesgourgues and S. Pastor, Phys. Rept. 429, 307 (2006) |arXiv:astro-ph/0603494|. 

[3] A. Strumia and F. Vissani, arXiv:hep-ph/0606054. 

[4] S. Hannestad, Ann. Rev. Nucl. Part. Sci. 56, 137 (2006) |arXiv:hep-ph/0602058|. 

[5] S. Eidelman et al [Particle Data Group], Phys. Lett. B 592, 1 (2004). 

[6] H. Abele, Prog. Part. Nucl. Phys. 60, 1 (2008). 

[7] R. N. Mohapatra, et al, Rep. Prog. Phys. 70, 1757 (2007). 



27 



[8] F. T. Avignone, S. R. Elliott and J. Engel, Rev. Mod. Phys. 80, 481 (2008) |a,rXiv:0708.TTm] 
[nucl-ex]]. 

[9] S. Weinberg, Cosmology, Oxford University Press, New York (2008). 
[10] S. Dodelson, Modern Cosmology, Academic Press, San Diego (2003). 
[11] P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75, 559 (2003) [arXiv:astro-ph/0207347j . 
[12] S. M. Carroll, eConf C0307282, TTH09 (2003) [AIP Conf. Proc. 743, 16 (2005)] 

[arXiv: astro-pli/03 10342 [ . 
[13] E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 (2006) 

[arX iv:hep-th/0603 057| . 
[14] S. Nojiri and S. D. Odintsov, Phys. Rev. D 68, 123512 (2003). 

[15] S. M. Carroll, A. De Fehce, V. Duvvuri, D. A. Easson, M. Trodden and M. S. Turner, Phys. 

Rev. D 71, 063513 (2005). 
[16] B. Ratra and P. J. E. Peebles, Phys. Rev. D 37, 3406 (1988) 
[17] C. Wetterich, Nucl. Phys. B 302, 668 (1988). 

[18] A. Linde, Lect. Notes Phys. 738, 1 (2008) [ arXiv:0705.0T64] [hep-th]]: Particle physics and 

Inflationary Cosmology, Chur, Switzerland: Harwood (1990) [arXiv: |hep-th / 0503203| . 
[19] S. Weinberg, Phys. Rev. D 77, 123541 (2008) [arXiv:0804.4291 [hep-th]]. 
[20] E. Mocchiutti et aZ., a?Xiv:0905.2551 [astro-ph.HE]. 
[21] E. A. Baltz et al, JCAP 0807, 013 (2008) [arXiv:0806.2911 [astro-ph]]. 

[22] K. Kohri, A. Mazumdar, N. Sahu and P. Stephens, Phys. Rev. D 80, 061302 (2009) 
[a^iv:0907.0622l [hep-phj] . 



[23] G. W. Anderson and S. M. Carroll, arXiv :astro-ph/9711288 
[24] M. B. Hoffman, arXiv:astro-ph/0307350. 

[25] G. R. Farrar and P. J. E. Peebles, Astrophys. J. 604, 1 (2004) [ arXiv:astro-ph/0307316j . 
[26] M. Kawasaki, H. Murayama and T. Yanagida, Mod. Phys. Lett. A 7, 563 (1992). 



[27] J. Garcia-Belhdo, Int. J. Mod. Phys. D 2, 85 (1993) [ arXiv:hep-ph/9205216 . 

[28] D. Comelh, M. Pietroni and A. Riotto, Phys. Lett. B 571, 115 (2003) [arXiv:hep-ph/0302080| . 

[29] U. Alam, V. Sahni, T. D. Saini and A. A. Starobinsky, Mon. Not. Roy. Astron. Soc. 354, 275 
(2004) [arX iv:astro-ph/03113 64|; L. Amendola, M. Gasperini and F. Piazza, JCAP 0409, 014 
(2004) |ar Xiv:astro-ph/0407573 l. M. Szydlowski, T. Stachowiak and R. Wojtak, Phys. Rev. D 
73, 063516 (2006) |arXiv:astro-ph705ir650| ; H. Li, B. Feng, J. Q. Xia and X. Zhang, Phys. 
Rev. D 73, 103503 (2006) |arXiv:astro-ph/0509272|. Z. K. Guo, N. Ohta and S. Tsujikawa, 
Phys. Rev. D 76, 023508 (2007) |arXiv:astro-ph/07020T5]. 



[30] U. Franca and R. Rosenfeld, Phys. Rev. D 69, 063517 (2004) [a rXiv:astro-ph/0308149| . 

[31] U. Franca, M. Lattanzi, J. Lesgourgues and S. Pastor, Phys. Rev. D 80, 083506 (2009) 
|^Xiv:0908.0534 [astro-ph.CO]. 

[32] G. Huey and B. D. Wandelt, Phys. Rev. D 74, 023519 (2006) [arXi v:astro-ph/0 407196|. 

[33] B. Wang, J. Zang, C. Y. Lin, E. Abdalla and S. Micheletti, Cosmological Parameters," Nucl. 
Phys. B 778, 69 (2007) [arXiv:astro-ph/0607126|. 

[34] R. Mainini and S. Bonometto, JCAP 0709, 017 (2007) [arX iv:0709.0174 l [astro-ph]]. 

[35] T. Koivisto, Phys. Rev. D 72, 043516 (2005) |arXiv:astro-ph /0504571|; R. Rosenfeld, Phys. 
Rev. D 75, 083509 (2007) | arXiv:astro-p h/0701213| ; S. Das, P. S. Corasaniti and J. Khoury, 
Phys. Rev. D 73, 083509 (2006) |arXiv:astro-ph/0510628]; X. J. Bi, B. Feng, H. Li and 
X. M. Zhang, neutrinos," Phys. Rev. D 72, 123523 (2005) [arXiv:hep-p h/04120 02|; M. Man- 
era and D. F. Mota, Mon. Not. Roy. Astron. Soc. 371, 1373 (2006) |ar Xiv:astro-p h/0504519|; 
L. Amendola, G. Camargo Campos and R. Rosenfeld, Phys. Rev. D 75, 083506 (2007) 



28 



|arXiv:astro-ph/0610806| ; G. J. Steph enson, J. T. Goldman and B. H. J. McKellar, Int. J. Mod. 
Phys. A 13, 2765 (1998) | arXiv:hep-ph /9603392;| ; S. Matarrese, M. Pietroni and C. Schimd, 
JCAP 0308, 005 (2003) |arXiv:astro-ph/0305224|; A. V. Maccio, C. Quercellini, R. Mainini, 
L. Amendola and S. A. Bonometto, Phys. Rev. D 69, 123516 (2004) | arXiv:astro-p h/0309671 ]; 



F. Vernizzi, Phys. Rev. D 69, 083526 (2004) |arXiv:astro-ph/0311167 



H. Li, M. Z. Li and 



X. M. Zhang, Phys. Rev. D 70, 047302 (2004) |arXiv:hep-ph/0403281|; K. Cheung and 



O. Seto, Phys. Rev. D 69, 113009 (2004) |arXiv:hep-ph/0403003 ; A. Nusser, S. S. Gubser 



and P. J. E. Peebles, Phys. Rev. D 71, 083505 (2005) |arXiv:astro-ph/0412586 1. S. S. Gubser 
and P. J. E. Peebles, Phys. Rev. D 70, 123510 (2004) |arXiv:hep-th/ 0402225| . 

[36] R. Fardon, A. E. Nelson and N. Weiner, JCAP 0410, 005 (2004) |a rXiv:astro-ph/0309800j . 

[37] P. Gu, X. Wang and X. Zhang, Phys. Rev. D 68, 087301 (2003). 



[38] R. D. Peccei, Phys. Rev. D 71, 023527 (2005) |arXiv:hep-ph/0411137 



[39] J. Grande, J. Sola and H. Stefancic, JCAP 0608, Oil (2006) [arXiv:gr-qc/0604057|; J. Grande, 
A. Pelinson and J. Sola, Phys. Rev. D 79, 043006 (2009) [arXiv:0809.3462 [astro-ph]]. 

[40] N. Afshordi, M. Zaldarriaga and K. Kohri, Phys. Rev. D 72, 065024 (2005) 
(arXiv : astro-ph /0 506663 1 . 

[41] M. Kaphnghat and A. Rajaraman, Phys. Rev. D 75, 103504 (2007) [arXiv:astro-ph/0601517|. 

[42] V. Pettorino and C. Baccigalupi, Phys. Rev. D 77, 103003 (2008) [arXiv: 0802. 1086 [astro-ph]]. 

[43] O. E. Bjaelde, A. W. Brookfield, C. van de Bruck, S. Hannestad, D. F. Mota, L. Schrempp 
and D. Tocchini-Valentini, JCAP 0801, 026 (2008) [arXi v:0705.2018 [astro-ph]]; O. E. Bjaelde 
and S. Hannestad, Phys. Rev. D 81, 063001 (2010) [ arXiv:()8()6.214 6 [astro-ph]]. 

[44] R. Bean, E. E. Flanagan and M. Trodden, Phys. Rev. D 78, 023009 (2008) [ar Xiv:()7()9.Tl28] 
[astro-ph]]; New J. Phys. 10, 033006 (2008) [arXiv:0709.1124 [astro-ph]]. 

[45] J. Valiviita, E. Majerotto and R. Maartens, JCAP 0807, 020 (2008) [a rXiv:0804.0232l [astro- 

Ph]]- 

[46] C. Wetterich and V. Pettorino, la,rXiv:0905.071 5 [astro-ph.CO]; V. Pettorino, D. F. Mota, 

G. Robbers and C. Wetterich, AIP Conf. Proc. 1115, 291 (2009) [arXiv:0901.1239 [astro-ph]]; 
L. Amendola, M. Baldi and C. Wetterich, Phys. Rev. D 78, 023015 (2008) [arXiv:0706.3064, 
[astro-ph]]; C. Wetterich, Phys. Lett. B 655, 201 (2007) [arXiv:0706.4427 [hep-ph]]; 
D. F. Mota, V. Pettorino, G. Robbers and C. Wetterich, Phys. Lett. B 663, 160 (2008) 
[arXiv:0802J5l5] [astro-ph]]; V. Pettorino, D. F. Mota, G. Robbers and C. Wetterich, AIP 
Conf. Proc. 1115, 291 (2009) [arXiv: 09()1.123 9 [astro-ph]]. 

[47] A. W. Brookfield, C. van de Bruck and L. M. H. Hall, Phys. Rev. D 77, 043006 (2008) 
[arXiv:0709.2297 [astro-ph]]; A. W. Brookfield, C. van de Bruck, D. F. Mota and D. Tocchini- 
Valentini, Phys. Rev. Lett. 96, 061301 (2006) |arXiv:astro-ph/0503349|; Phys. Rev. D 73, 
083515 (2006) [Erratum-ibid. D 76, 049901 (2007)] |arXiv:astro-ph705 12367] . 

[48] A. E. Bernardini and O. Bertolami, Phys. Rev. D 77, 083506 (2008) [arXiv:0712.1534 [astro- 
ph]]; Phys. Lett. B 662, 97 (2008) [arXiv: 0802. 4449 [hep-ph]]; Phys. Lett. B 684, 96 (2010) 
[arXiv:0909.1280 [gr-qc]]; Phys. Rev. D 80, 123011 (2009) [arXiv:090 9.1541 [gr-qc]]. 

[49] R. Takahashi and M. Tanimoto, Phys. Lett. B 633, 675 (2006) [arXiv:hep-ph/0507142] . 



[50] R. Takahashi and M. Tanimoto, JHEP 0605, 021 (2006) |arXiv:astro-ph/0601119 



[51] R. Fardon, A. E. Nelson and N. Weiner, JHEP 0603, 042 (2006) |arXiv :hep-ph/ 0507235l. 
[52] C. Spitzer, arXiv:astro-ph/0606034'. 



[53] K. Ichiki and Y. Y. Keum, JCAP 0806, 005 (2008) [a rXiv:0705.213 4 [astro-ph]]; JHEP 0806, 

058 (2008) [arXiv:0803.2274 [hep-ph]]. 
[54] J.I. Kapusta and C. Gale, Finite- Temperature Field Theory: Principles and Applications, 



29 



Second Edition, Cambridge University Press, Cambridge (2006). 
[55] J. Zinn- Justin, Quantum Field Theory and Critical Phenomena, 4th edition. Clarendon Press, 
Oxford (2002). 

[56] W. Greiner, L. Neise, and H. Stocker, Thermodynamics and Statistical Mechanics, Springer- 
Verlag, New York (1995). 



[57] S. Hannestad, JCAP 0305, 004 (2003) [ arXiv:astro -ph70303076 ] ; V. Barger, J. P. Kneller, 
H. S. Lee, D. Marfatia and G. Steigman, Phys. Lett. B 566, 8 (2003) [arXiv :hep-ph/0305075| 
A. Melchiorri and C. J. Odman, Phys. Rev. D 67, 081302 (2003) |arXiv:astro-ph/0302361 ] 
S. D. Stirhng and R. J. Scherrer, Phys. Rev. D 66, 043531 (2002) |arXi v:astro-ph/020 6173| 
J. Hamann, J. Lesgourgues and G. Mangano, JCAP 0803, 004 (2008) [a rXiv:0712.2826l [astro- 

Ph]]. 

[58] L L. Shapiro and J. Sola, Phys. Lett. B 475, 236 (2000) [arXiv:hep-ph/9910462|; JHEP 0202, 
006 (2002) [arXiv:hep-th/0012227l; L L. Shapiro, J. Sola, C. Espana-Bonet and P. Ruiz- 
Lapuente, Phys. Lett. B 574, 149 (2003) [arXiv: astro-ph/0303306 ] ; I. L. Shapiro, J. Sola and 
H. Stefancic, JCAP 0501, 012 (2005) |arXiv:hep-ph/0410095| V J. Sola, J. Phys. A 41, 164066 
(2008) (a,rXiv:07l 0.4151 1 [hen-th]]. 
[59] F. Bauer, J. Sola and H. Stefancic, Phys. Lett. B 678, 427 (2009) iarXiv:0902.2215 [hep-th]]; 
Phys. Lett. B 688, 269 (2010) [arXiv:0912.0677 [hep-th]]; D. A. Demir, Found. Phys. 39, 
1407 (2009) [arXiv:0910.2730 [hep-th]]. S. Basilakos, M. Plionis and J. Sola, Phys. Rev. D 80, 
083511 (2009) [arXiv:0907.4555l[astro-ph.CO]]. 
[60] N.D. Birrel and P.C.W. Davies, Quantum Fields in Curved Space, Cambridge University Press, 
Cambridge (1982). 

[61] G. E. Volovik, Phil. Trans. Roy. Soc. Lond. A 366, 2935 (2008) [arXiv:0801.0724 [gr-qc]]; 

F. R. Khnkhamer and G. E. Volovik, Phys. Rev. D 77, 085015 (2008) [arXiv:0711.3170 [gr- 
qc]]; G. E. Volovik, The Universe in a Helium Droplet, Clarendon Press, Oxford (2003); 

G. E. Volovik, JETP Lett. 80, 531 (2004); JETP Lett. 77, 769 (2003). 
[62] H. J. de Vega and N. G. Sanchez, arXiv:astro-ph/0701212| 

[63] E. V. Gorbar and L L. Shapiro, JHEP 0302, 021 (2003) [arXiv:hep-ph/0210388|; JHEP 0306, 
004 (2003) [arXiv:hep-ph/030312 4|; L L. Shapiro and J. Sola, Phys. Lett. B682, 105 (2009) 
|arXiv:091 0.4935^hep-th]] . 

[64] A detailed analysis of the vicinity of the critical point shows that it is an unstable degenerate 
saddle-node of the differential equation ()107p . See, e.g., N.N. Bautin and E.A. Leontovich, 
Methods of the Qualitative Analysis of Dynamical Systems on a Plane, Second Edition, Nauka, 
Moscow (1990). 

[65] E. E. O. Ishida, R. R. R. Reis, A. V. Toribio and L Waga, Astropart. Phys. 28, 547 (2008) 

|arXiv:0706.0546] [astro-ph]]. 
[66] S. M. R,. Micheletti. larXiv:1009.6198l [gr-ac]. 



